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The infrared properties of lattice Landau gauge QCD of SU(3) are studied by measuring gluon propagator, ghost 
propagator, QCD running coupling and Kugo-Ojima parameter of /3 = 6.0, le"*, 24'', 32" and = 6.4, 32'', 48'', Se"* 
lattices. By the larger lattice measurements, we observe that the runnning coupling measured by the product 
of the gluon dressing function and the ghost dressing function squared rescaled to the perturbative QCD results 
near the highest lattice momentum has the maximum of about 2.2 at around q = 0.5 GeV/c, and behaves either 
approaching constant or even decreasing as q approaches zero. The magnitude of the Kugo-Ojima parameter 
is getting larger but staying around —0.83 in contrast to the expected value —1 in the continuum theory. We 
observe, however, there is an exceptional sample which has larger magnitude of the Kugo-Ojima parameter and 
stronger infrared singularity of the ghost propagator. The reflection positivity of the 1-d Fourier transform of the 
gluon propagator of the exceptional sample is manifestly violated. 

Gribov noise problem was studied by performing the fundamental modular gauge (FMG) fixing with use of the 
parallel tempering method of /3 = 2.2, 16 SU(2) configurations. Findings are that the gluon propagator almost 
does not suffer noises, but the Kugo-Ojima parameter and the ghost propagator in the FMG becomes ~ 5% less 
in the infrared region than those suff'ering noises. It is expected that these qualitative aspects seen in SU(2) will 
refiect in the infrared properties of SU(3) QCD as well. 



1. Introduction 

One of our basic motivations in the present 
study is verification of the color confinement 
mechanism in the Landau gauge. Two decades 
ago, Kugo and Ojima proposed a criterion for 
the color confinement in Landau gauge QCD 
using the Becchi-Rouet-Stora-Tyutin(BRST) in- 
variance of continuum theory pP . Gribov pointed 
out that the Landau gauge can not be uniquely 
fixed, that is, the Gribov copy problem, and ar- 
gued that the unique choice of the gauge copy 
could be a cause of the color confinement [2]. 
Later Zwanziger developed extensively the lat- 
tice Landau gauge formulation (51 in view of the 
Gribov copy problem. Kugo and Ojima started 
from naive Faddeev-Popov Lagrangian obviously 
ignoring the Gribov copy problem, and gave the 
color confinement criterion with use of the foUow- 
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ing two point function of the covariant derivative 
of the ghost and the commutator of the antighost 
and gauge field: 

x,y \ ^ / xy 

where lattice simulation counterpart is utilized. 
They claim that sufficient condition of the color 
confinement is that m(0) = —1 with u"-^{Q) = 
6''''u{<d). Kugo showed that 

l + .(0) = |i = |, (2) 

where Z^ is the gluon wave function renormaliza- 
tion factor, Zi is the gluon vertex renormalization 
factor, and Z^ is the ghost wave function renor- 
malization factor, respectively. In the continuum 
theory Zi is a constant in perturbation theory 
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and is set to be 1. On the lattice, it is not evident 
that it remains 1 when strong non-perturbative 
effects are present. In a recent SU(2) lattice sim- 
ulation with several values of (3, finiteness of Zi 
seems to be confirmed, but its value differ from 
1. 

The same equality, the first one of equations 
(O, was derived by Zwanziger with use of his 
"horizon condition" about the same time ,3_. It 
is to be noted that arguments of both Kugo and 
Zwanziger are perturbative ones in that they used 
diagramatic expansion, and the equation (O is of 
continuum theory or continuum limit. 

The non-perturbative color confinement mech- 
anism was studied with the Dyson-Schwinger ap- 
proach [4.5 and lattice simulations l 6l7l8l9lfU| . 
Both types of studies are complementary in that 
Dyson-Schwinger approach needs ansatz for trun- 
cation of interaction kernels and lattice simula- 
tion is hard to draw conclusions of continuum 
limit although the calculation is one from the first 
principle. 

We measured in SU(3) lattice Landau gauge 
simulation with use of two options of gauge 
field definition (log?7, U linear; see below), 
gluon propagator, ghost propagator, QCD run- 
ning coupling and Kugo-Ojima parameter of 
13 = 6.0,16^24^324 and (3 = 6.4, 32^, 48^ 56^ 
lattices. The QCD running coupling Us = 
(7^/47r can be measured in terms of gluon dress- 
ing functiuon Z{q^) and ghost dressing function 
G{q'^), as renormalization group invariant quatity 
g'^G{q^Y Z{q^). Infrared features of is not 
known, however, and there remains a problem 
of checking the Gribov noise effect among those 
quantities, since there exist no practical algo- 
rithms available so far for fixing fundamaental 
modular gauge (see below). In our 56"* simula- 
tion, we encountered a copy of an exceptional 
configuration yielding extraordinarily large Kugo- 
Ojima marameter c — — u(0), and studied its fea- 
ture in some more detail. We made another copy 
by adjusting controlling parameter in gauge fix- 
ing algorithm, and measured copywise 1-d FT of 
the gluon propagator, and found violation of re- 
flection positivity in both cases. 

In order to study the Gribov copy problem, we 
made use of parallel tempering 9 with 24 replicas 



to fix fundamental modular gauge in SU(2), [3 = 
2.2, 16^ lattice, and obtained qualitatively similar 
result as Cucchieri)12|. 

1.1. The lattice Landau gauge 

We adopt two types of the gauge field defini- 
tions: 

1. logU type: U^^^ = e^--", ^ = -A^,^,, 

2. U linear type: ^j;,^ l^{Ux,^l - Ul ^^)\trip., 

where trip, implies traceless part. The Landau 
gauge, dA^ = is specified as a stationary point 
of some optimizing functions Fu{g) along gauge 
orbit where g denotes gauge transformation, i.e., 
6Fu{g) = for any 6g. 

Here Fu{g) for the two options are |11)I3) 

1. Fuig) = WA^r = E.,^tr [a^^JaQ, 

2- Fuig) = E.,,. (1 - trUiJ , 

respectively. Under infinitesimal gauge transfor- 
mation g~^Sg = e, its variation reads for either 
defintion as 

AFuig) = -2{dAS\e) + (e| - dD{US)\e) + ■■■, 

where the covariant derivativativc D^{U) for two 
options reads commonly as 

where df^cf) — (jj^x + fi) — (j){x), and 

7 0(a; + /i) +(t){x) , . , n ■.■ r 

(p — , but definition of operation 

S{Ux,^)B.j;^i_, is given 

1. S{Ux,fj.)Bx.fj. = T{Ax.ji)Bx^^ 
where Ax,i_i = adjA^^^ = [Ax,ii, ■], 



(3) 



T(x) 



x/2 
th(a;/2)' 



2. S{Ux,fj,)Bx^fj, — — 



1 Ux,^ + 



B 



(4) 



trip. 

The fundamental modular gauge (FMG) 
[31 is specified by the global minimum of the 

minimizing function Fu{g) along the gauge or- 
bits in either case, i.e., 

A ^ {U\ A ^ AiU),Fu{l) - MuigFuia)}, 
A C ft, where is Gribov region (local min- 
ima), and n = {U\- dD{U) > , dA^O}. 
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2. Numerical simulation of lattice Landau 
gauge QCD 

2.1. Method of simulation 

First we produce Monte Carlo (Boltzmann) 
samples of link variable configuration according 
to Wilson's plaquette action by using the heat- 
bath method. We use optimally tuned combined 
algorithm of Creutz's and Kennedy-Pendleton's 
for SU(2) case, and Cabbibo -Marinari pseudo- 
heat-bath method with use of the above SU(2) 
algorithm for SU(3) case. As for FMG fixing, 
the smearing gauge fixing works well for large (3 
on relatively small size lattice 8 . However we 
found that its performance for SU(3), (3 — 6.0, 
16'* lattice is not perfect in comparison with our 
standard method for log U type definition, in 
which the third order perturbative treatment of 
the linear equation with respect to gauge field A, 
^dDe = dA is performed. For U linear defini- 
tion, we use the standard overrelaxation method 
for both SU(2) and SU(3). Thus we know that 
our gauge fixing is not FMG fixng. The accuracy 
of dA{U) = is 10^^ in the maximum norm in 
all cases. 

Measurement of quantity in question, gluon 
propagator, ghost propagator, Kugo-Ojima pa- 
rameter, etc., is performed with use of Landau 
gauge fixed copy, and averaged over Monte Carlo 
samples. 

We define the gluon dressing function Z^(g^) 
from the gluon propagator of SU(n) 

DM = E e-''=^tr(A^(x)t^.(0)) 

Tl i 

= iS,^-^^)DMq'), (5) 
as ZA{q^) = q^DAiq^). 

The ghost propagator is defined as the Fourier 
transform of an expectation value of the inverse 
of Faddeev-Popov(FP) operator A4 = —dD, 

D^a\^,y) = {tr(\\\{M[U])-'\\''y)), (6) 

where the outmost (• ■ •) denotes average over 
samples U. 

2.2. Numerical results of simulation 

DA{q^) is defined in 0, and its dressing func- 
tion is given as Za((7^) = q^DA{q^)- Numerical 
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Figure 1. The gluon dressing function as the 
function of the momentum g(GeV). /3 = 6.4, 
48* (stars) and 56** (diamonds) in the logU defini- 
tion. The solid line is that of the MOM scheme. 



results are given in figure^ 

The ghost propagator is defined in ®. It 
should be noted that the operator [A^(t/)]!J^ 
is a real symmetric operator when dA ~ 0, 
and it transforms covariantly under global color 
gauge transformation g, i.e., [Ai{U^)]'^'' = 
[g^ {M{U)}g]'^''. The conjugate gradient method 
is used for the evaluation the ghost propagator, 
and absence of color off-diagonal components 
manifests itself, a signal of no color summetry 
violation in Landau gauge. The accuracy of the 
inversion is maintained less than 5% in maximum 
norm in error check of the source term. 

We measured the running coupling from the 
product of the gluon dressing function and the 
ghost dressing function squared. 

as{q') = ^Z{q')G{q'f ^ (ga)-^("°+2"-). (7) 

The lattice data (diamonds) are compared with 
a fit of the DSE approach 5 with infrared fixed 
point OfQ = 2.5, and the pertubative QCD with 
c/ q^ correction [7] and the contour improved per- 
turbation methodfTU]. The The result using the 
ghost propagator of I a copy (star) is also plotted 
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Figure 2. The running coupling as{q) as a func- 
tion of the logarithm of momentum logj^Q g(GeV) 
of the log U, P = 6.4, 56^ lattice using the ghost 
propagator of the average (diamond). The result 
using the ghost propagator of Ia copy (star) is 
also plotted for comparison. The DSE approach 
(long dashed line), the perturbative QCD+c/g^ 
(short dashed line) and the contour improved per- 
turbation method (dotted line) are also shown. 



for presenting the ambiguity due to the Gribov 
copy. 

Kugo-Ojima color confinement criterion is 
given with use of u°'''{p'^) defined in ^ with 
M'»f'(p2) ^ S''''u{p^), as u(0) ^ -c- c ^ 1 ^ 
color confinement. Zwanziger's horizon condition 
[3] derived in the infinite volume limit is written 
as follows; 




where e = (E.,^ tr(A«t^([/^_^)A-)^ /{(n^ - 

1)1^}, and the horizon condition reads 
lim Gnu (p) — e = 0, and the l.h.s. of the con- 

dition is ^— ^ + {d — l)c ~ e ~ (d ~ l)h where 

e . . 

h = c and dimension d = 4, and it follows that 

d 



h = —>■ horizon condition, and thus the horizon 
condition coincides with Kugo-Ojima criterion 
provided the covariant derivative approaches the 
naive continuum limit, i.e., e/d = 1. 

Values of c, —,h are shown in the Table ^ (sufhx 
d 

1: logC/, suffix 2: L/— linear) 

Among gauge fixed configurations, /3 — 6.4, 
\ogU type, there appeared a configuration (Ia) 
which have exceptionally large Kugo-Ojima pa- 
rameters, c = 0.851(77) with rather large fiuc- 
tuation among diagonal components. For more 
informtion, we made another copy of it by chang- 
ing some parameter of the gauge fixing program. 
They have very close values in norm, but 

fairly large difference in physically important val- 
ues. One dimensional Fourier transforms trans- 
verse to 4 axes were calculated with respect to 
sample-wise transverse gluon propagator func- 
tions of momentum q. Existence of negative value 
region in one dimensional Fourier tranform of 
(sample averaged) propagator implies violation of 
refiection positivity postulate 

3. The Gribov copy problem 

The Gribov copy problem is one of fundamen- 
tal problems in QCD. FMG is the unique gauge 
known as the legitimate Landau gauge without 
Gribov copies. No algorithms for FMG fixing 
have been established. Numerical invesigations of 
Gribov noise have been done in a few works jl2| . 
For U linear type gauge field definition, FMG 
gauge fixing with use of parallel tempering (FT) 
with 24 replicas was performed for /? = 2.2, 16''', 
SU(2) 67 Monte Carlo configurations. Some more 
details of the method are given in |5] . We com- 
pared data between two groups, the absolute min- 
imum configurations of Fij{g) obtained by Lan- 
dau gauge fixing via PT and the 1st copies ob- 
tained by the direct application of Landau gauge 
fixing on Monte Carlo configurations. The re- 
sults are as follows; l)The singularity of the ghost 
propagator of the FMG is about 6% less than that 
of the 1st copy. 2)Kugo-Ojima parameter c of the 
FMG is about 4% smaller than that of the 1st 
copy. 3) The gluon propagator of the two groups 
are almost the same within statistical errors. 4) 
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Table 1 



The Kugo-Ojima parameter c in logU and J7— linear version. (3 — 6.0 and 6.4. 



f3 


L 


Cl 


ei/d 


hi 


C2 


e2/d 


h2 


6.0 


16 


0.628(94) 


0.943(1) 


-0.32(9) 


0.576(79) 


0.860(1) 


-0.28(8) 


6.0 


24 


0.774(76) 


0.944(1) 


-0.17(8) 


0.695(63) 


0.861(1) 


-0.17(6) 


6.0 


32 


0.777(46) 


0.944(1) 


-0.16(5) 


0.706(39) 


0.862(1) 


-0.15(4) 


6.4 


32 


0.700(42) 


0.953(1) 


-0.25(4) 


0.650(39) 


0.883(1) 


-0.23(4) 


6.4 


48 


0.793(61) 


0.954(1) 


-0.16(6) 


0.739(65) 


0.884(1) 


-0.15(7) 


6.4 


56 


0.827(27) 


0.954(1) 


-0.12(3) 


0.758(52) 


0.884(1) 


-0.13(5) 



The running coupling of FMG is about 10% sup- 
pressed than that of the 1st copy. 5) The horizon 
function deviation parameter h of the FMG is not 
closer to 0, i.e. the value expected in the contin- 
uum limit, than that of the 1st copy. 

In Figl^l shown is the scatter plot of Gribov 
copies on the plane of Kugo-Ojima parameter vs 
1 — Fu(g), 32 copies via random gauge trans- 
formation and one copy given via FT, all being 
Gribov copies on a single gauge orbit. 
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Figure 3. The Kugo-Ojima parameter w(0) vs 
1 — Fij{g) of the 33 copies on a single gauge or- 
bit. S'[/(2),@/3 = 2.2, 16". The maximum opti- 
mization point by Cucchieri's prescription(lower 
circle) and our prescription(upper circle) are in- 
dicated. 



Figure 4. The Kugo-Ojima parameter u(0) vs 
1 — Fu{g) of 16 Gribov copies via FT for three 
gauge orbits each. 



In Fig^ the best one of 16 copies should be 
picked up for the FMG. This figure shows how 
flat is the bottom of Fu{g) with respect to the 
physical quantity. 

4. Discussion and Conclusion 

The SU(3) 13 = 6.0, 16^ 24", 32" and /? = 
6.4, 32'*, 48'', 56" lattice Landau gauge simulation 
was performed. A comparison of the data of log U 
type and the t/— linear type reveals that the gluon 
propagator does not depend on the type and the 
ghost propagator of the ?7— linear is about 14% 
larger than the log U . The Kugo-Ojima parame- 
ter c is still getting larger w(0) ~ —0.83 in logf 
and about -0.76 in [/—linear, in accordance to 
lattice size. The running coupling as{q) of the 
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log U type has a maximum of about 2.2, and that 
of the J7— hnear is sHghtly smaller since they are 
rescaled at the high lattice momentum region. 
The absolute value of the index ac increases as 
the lattice size becomes large but since it is de- 
fined at q ~ OAGeV region, it is smaller than 
K ~ 0.5 of the DSE which is defined at the 
momentum region by about factor 2. 

As is studied in SU(2) FMG data, the gluon 
propagator suffers almost no Gribov noise, but 
magnitude of Kugo-Ojima parameter becomes 
smaller in FMG, and the ghost propagator ex- 
hibits less singular infrared behavior in FMG than 
that supposed suffering noise. Accordingly the 
running coupling becomes smaller in FMG. All 
samples in SU(3) can not be supposed in FMG 
and supposed to suffer the Gribov noise, and it 
is expected that these qualitative aspects seen 
in SU(2) will reflect in the infrared behavior of 
SU(3) QCD as well. The FMG is mathemat- 
ically well defined on lattice, and its existence 
can be proven. But rather flat valued feature 
of the FMG optimizing function with respect to 
physical quantities seem to keep annoying us, at 
least, technically, and this fact may suggest ne- 
cessity of some new formulation of infrared dy- 
namics of QCD. In the Langevin formulation of 
Landau gauge QCD, Zwanziger conjectured that 
the path integral over the FM region will become 
equivalent to that over the Gribov region in the 
continuum If his conjecture is correct, then 
in order for the gauge dependent quantity, e.g., 
the optimizing function to have the same expecta- 
tion values on both regions, the situation that the 
probability density may be concentrated in some 
common local region becomes favorable one. The 
proximity of the FM region and the boundary of 
the Gribov region in SU(2) in 8*, 12^ and 16^ lat- 
tices with [3 — 0,0.8,1.6 and 2.7 was studied in 
|14| . The tendency that the smallest eigenvalue 
of the Faddeev-Popov matrix of the FMG and 
that of the 1st copy come closer as /3 and lattice 
size become larger was observed, although as re- 
marked in 55 the physical volume oi (3 = 2.7, 
16^ lattice is small and not close to the contin- 
uum limit. In this respect, further study of lattice 
size and f3 dependence of Gribov noise of various 
quantities will become an important problem in 



future. 

Our observation of exceptional conflgurations 
in SU(3), f3 — 6.4, 56'* may be considered as a 
signal of the tendency of the probabililty concen- 
tration to the Gribov boundary 3 11. 

The confinement scenario was recently re- 
viewed in the framework of the renormalization 
group equation and dispersion relation !1^. It 
is shown that the gluon dressing function satis- 
fies the superconvergence relation, and the gluon 
propagator does not necessarily vanish as Gri- 
bov and Zwanziger conjectured, but it should be 
finite. The multiplicative renormalizable DSE 
approach [H] predicts that the exponent k ~ 0.5 
and the infrared fixed point — 2.6 and our 
lattice data are consistent with this prediction. 
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